On the detectability of quantum radiation in Bose-Einstein condensates 
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Based on doubly detuned Raman transitions between (meta) stable atomic or molecular states 
and recently developed atom counting techniques, a detection scheme for sound waves in dilute 
Bose-Einstein condensates is proposed whose accuracy might reach down to the level of a few 
or even single phonons. This scheme could open up a new range of applications including the 
experimental observation of quantum radiation phenomena such as the Hawking effect in sonic 
black-hole analogues or the acoustic analogue of cosmological particle creation. 
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PACS numbers: 04.70.Dy, 03.75.Kk, 42.65.Dr. 

Gaseous atomic or molecular Bose-Einstein conden- 
sates are in several ways superior to other superfLuids: 
Apart from a very good theoretical understanding via the 
mean- field formalism (in the dilute-gas limit), they of- 
fer unprecedented options for experimental manipulation 
and control. It is possible to influence the shape, den- 
sity, flow profile, and coupling strength of Bose-Einstein 
condensates via external electromagnetic fields. Finally, 
these condensed gases are rather robust against the im- 
pact of the environment such that one may reach ex- 
tremely low temperatures. 

In view of all these advantages, the question naturally 
arises whether it could be possible to measure so far un- 
observed quantum radiation phenomena in a suitable set- 
up. These exotic quantum effects include cosmological 
particle creation (due to the amplification of quantum 
fiuctuations in an expanding/contracting universe 
as well as the acoustic analogue of Hawking radiation 
HQ in "dumb holes" [IIIH- 

For wavelengths which are much longer than the heal- 
ing length ^, the propagation of phonons in Bose-Einstein 
condensates is analogous to a scalar field in a curved 
space-time described by the effective metric 
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which is determined by the density qq and velocity Vq of 
the background fluid. For example, assuming an effec- 
tively one-dimensional stationary flow, the point where 
the fluid velocity exceeds the local speed of sound Cg 
corresponds to the sonic analogue of the horizon of a 
black hole. The corresponding Hawking temperature is 
determined by the velocity gradient 
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i.e., the characteristic length scale A over which the flow 
changes. Since this length scale should be large compared 
to the healing length ^ (typically of order micrometer) 
for the curved space-time analogy to apply, a speed of 
sound of order mm/s leads to an upper bound for the 
typical energy of the Hawking phonons of order 10~^^ eV 
corresponding to a temperature on the nano-Kelvin level. 



Moreover, since the fluid velocity equals the sound 
speed at the acoustic horizon, only a limited number of 
these low-energy phonons will be created by the Hawking 
effect - unless one has a very large reservoir for the con- 
densate flow: Since the Hawking radiation is thermal, the 
typical distance between two emitted Hawking phonons 
is given by their characteristic wavelength A and hence 
it is much larger than the healing length ^. In addi- 
tion, Bose-Einstein condensates are formed by atoms (or 
molecules) whose inter-particle distance ad is far bigger 
than the s-wave scattering length (dilute-gas limit). 
As a result, a healing length ^ oc a^-y/ ad/cLs contains 
many atoms ^ S> ad, i-e., we have a hierarchy of length 
scales A 3> ^ S> fld 3> fls. Consequently, the number 
of Hawking phonons is extremely small compared to the 
number of atoms in the condensate a^/A 1. 

Similar arguments apply to the analogue of cosmologi- 
cal particle creation, which require a non-stationary set- 
up. Considering the effective metric in Eq. there are 
basically two possibilities for simulating the cosmic ex- 
pansion in Bose-Einstein condensates: an expansion of 
the condensate or a temporal variation of the speed of 
sound (which can be achieved via varying a^ by means of 
a Feshbach resonance, for example). For simplicity, we 
shall focus on the second possibility in the following, but 
the general ideas apply to both scenarios. The typical 
wavelength A of the created phonons is determined by 
the rate of change A = 0{c^g/cs) of Cs and should again 
be large compared to the healing length ^ for the curved 
space-time analogy to apply. In the absence of ampli- 
fication mechanisms such as resonances, the number of 
created phonon per wavelength is again of order one 0|. 

For a small number of phonons with an energy of or- 
der 10"^'^ eV, the usual detection mechanisms for sound 
via mechanical motion are extremely difficult to apply 
since the kinetic energy of a single atom with a veloc- 
ity of order mm/s already yields this amount. Usually, 
these measurements involve many atoms and thus many 
phonons (limit of classical waves). For example, it was 
possible to excite phonon modes via light scattering 
and to map out the dispersion relation etc. An indirect 
observation of the phonon number was achieved in ultra- 
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sensitive temperature measurements joj, which reached 
very low energies. 

Fortunately, a non-mechanical detection mechanism 
may circumvent these obstacles. For example, the ion- 
trap quantum computer in ^3 is based on optically in- 
duced transitions which require the simultaneous absorp- 
tion of a phonon in a given mode (due to the detuning 
of the Laser). Using the occurrence of the transition and 
the photon emitted during the decay back to the ground 
state as an indicator for the existence of the phonon yields 
an energy amplification over many orders of magnitude; 
and the detection of single photons in the optical range is 
difficult but feasible (in principle). A further amplifica- 
tion is possible if the phonon-assisted transition mediates 
between (meta) stable atomic states which can be sepa- 
rated or addressed individually: The detection of a small 
number of atoms can be achieved via fluorescence mea- 
surements involving many photons |llL Il2| . 
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FIG. 1: Sketch (not to scale) of the three-level (A) system 
and the doubly detuned Raman transitions denoted by Qi^2- 

In the following, a scheme for the transformation of 
low-energy phonons in a given mode into an equal num- 
ber of atoms in a different atomic state with controlled 
energy and momentum based on doubly detuned Raman 
transitions is presented. Let us consider atoms which can 
be described by a three- level (A) system consisting of two 
(meta) stable states and together with a third ex- 
cited level 5*3 with the energies cui < uj2 < ^^3- This 
three-level system is illuminated by two optical Laser 
beams which consist of many photons and can there- 
fore be treated as rapidly oscillating classical fields de- 
scribed by the effective Rabi frequencies fti{t) and 5^2 (i). 
Within the rotating wave and dipole approximation, the 
Lagrangian reads {h = I) 
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i'^'li'l + i*2*2 + - - t^2|^'2l 

-Losi^ll + [niit)^i^3 + n2it)%^3 + H.c] . (3) 



The frequencies of the two doubly detuned Laser beams 
are chosen according to (see Fig. ^ 
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with a large detuning A and a small detuning S (which 
will later determine the phonon energy). Introduc- 
ing the slowly varying variables ipi, ip2, and ips via 



*i(<) = Vi(i)exp{-iwi<}, *2(i) = V'2(i) exp{-iw20: 
and 'ipiit) exp{— i(a;3 -|- A)t}, we may solve the equation 
for the upper level V'a approximately for large detuning A 
via ^3 = -{^1^1 + nie-^^^2)/A + 0(1/A2). Insertion 
into Eq. Q yields the Lagrangian for the remaining two 
levels in the adiabatic approximation \il)3\ <C A 
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Assuming = \il2\ = il, both levels acquire the same 
additional shift fl^/A; otherwise we would obtain an ef- 
fective detuning S ~> 6' shifted by (|f7?| - \fll\)/A. 

An ideal quantum gas containing many of these atoms 
with mass m can be described by the many-particle field 
operator with the dynamics (Heisenberg picture) 
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where r,s = 1,2 are labels for the remaining two levels 
and Vr the corresponding potentials. The anti-hermitian 
space-time dependent transition amplitude Si2(t, x) = 
exp{—iSt + iK ■ r}fl^/ A represents the mode-coupling in 
Eq. where k arises from a small angle between the 
Raman beams and the resulting wavenumber mismatch 

_ .Laser .Laser 
fvi — tti-^ "''2, 

An expansion into single-particle energy-eigenstates 
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diagonalizes Eq. ® apart from the transitions, which are 
(in the rotating wave approximation) only relevant for 
Ei^a ~ ^2, ,3 = 6 (energy conservation) and if the spatial 
matrix element (/ia| S12 I/2/3) is large enough. For nearly 
homogeneous potentials Vr ~ const, the eigenfunctions 
are plane waves a k with Er.k = fc^/(2m) + Vr and 
the latter condition represents momentum conservation 
K = ki — k2. Hence, for a given frequency and wavenum- 
ber mismatch of the Lasers (i5, k), these energy and mo- 
mentum conservation conditions determine fci and fe2 up 
to a contribution perpendicular to n. For effectively one- 
dimensional condensates, therefore, we can address single 
modes ki = kie^ and k2 — ^263; by adjusting the Lasers. 

Now let us consider the following gedanken experi- 
ment: Initially all atoms are in the state r = 1 and form 
a nearly homogeneous and (quasi) one-dimensional con- 
densate, which is not in its ground state but contains a 
single phonon with a given wavenumber fcp — kpe^. In 
contrast to Eq. ((HJ, this requires a non- vanishing cou- 
pling g. However, if we switch off this interaction g 
adiabatically (e.g., via a Feshbach resonance), the sys- 
tem stays in this first excited state and finally contains 
a single atom with the momentum fcp of the original 
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phonon. After applying a Raman 7r-pulse (with the du- 
ration T ~ 7rA/57^) adapted to this wavenumber, e.g., 
K = kp and S = kp/{2m) + Vi — V2, exactly this single 
atom will be transferred to the other state r — 2, while 
all the condensate atoms are not affected (assuming that 
rotating wave approximation applies). 

If we can separate the two species r — 1 and r = 2 or 
address them individually, the number of atoms in the 
state r = 2 can be counted via fluorescence measure- 
ments [H Ei and yields ( in the ideal case) the number 
of phonons in a given mode kp present initially, i.e., one. 
For example, a beam with a frequency just between the 
resonances of the two species r = 1 and r = 2 is repulsive 
for one component and attractive for the other one and 
could be used as optical tweezers. 

With fixed momentum kp (homogeneous condensate), 
the energy gap between the ground state and the one- 
particle excited state decreases with diminishing cou- 
pling g in view of the Bogoliubov dispersion relation 
Lu'^{k) ~ ggk^/m + fc^/(2m)^. Hence, let us study the 
application of the Raman transitions in the presence of a 
non-vanishing coupling g (respecting the altered disper- 
sion relation). With interactions, the field operator i/^r of 
the two levels r = 1 , 2 obeys the equation of motion 
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For simplicity, we assume equal one-particle trapping po- 
tentials Vi —V2 = V and coupling constants 511 = (722 = 
<7i2 = 521 = g for the two species (otherwise we would 
again obtain an effective detuning S S'). Initially, all 
the atoms (and hence also the condensate) are in the 
state r — I, which facilitates the mean- field expansion 
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with the condensate wave-function ipc and the one- 
particle excitations x and The operator N — A'A 
counts the total number of particles. 

In complete analogy to the previous example, the con- 
densate ipc is not affected by the Raman transitions for 
S > (assuming that the rotating wave approximation 
applies) and hence the dynamics of the excitations read 
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with S — S12 — exp{—i(5i+iK-r}f2^/A as before. Assum- 
ing a nearly homogeneous condensate with + 5 ["iAcI ~ 
fi = const, a normal mode expansion yields 
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The atomic one-particle excitation operator Xk can be 
decomposed into phonon creation and annihilation oper- 
ators a\, and a^, respectively (m — 1 for brevity) 
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Inserting the time-dependences Cfe(^) = Cfce^''^^'" 7(2™)]* 
and ak{t) = afce~'"fc* as well as applying the rotating 
wave approximation, only the second term survives 




and we obtain the expected resonance condition 



(5 = cjfc 



(fc + Kf 

2m 



(13) 



(14) 



which implies energy conservation. Of course, for 
K?,k^ ^ 1/C^7 we reproduce the previous result 
Far below the healing length K^,k^ ^ 1/^'^ — m,gg, 
i.e., in the phonon regime, Eq. (|14|l simplifies to S — LUk 
and we get idCk+n/dt = fi/6ak f2^/A. Note that 
the pre-factor \/Jjj5 is a consequence of the interactions 
and illustrates the difference between phonons ( "dressed" 
atoms) and free-particle excitations as in Eq. 

The effective interaction Hamiltonian from Eq. (|13|) . 



int 



2ujk \2 fcV 



(cl 



,afc + CfcH 



,(15) 



has the following intuitive interpretation: Due to the de- 
tuning of the Raman beams, the missing energy 5 pro- 
hibits transitions from the multi-particle ground state of 
the condensate (which has a sharp and well-defined en- 
ergy) in component r = 1 to the other level r = 2 and 
must be compensated by absorbing a phonon with this 
(for k = —k) or an even higher energy ujk > 

If there are n phonons to annihilate (ofc), n atoms can 
be transferred to the r = 2 state (Cfc+K,) such that the 
final number of these transferred atoms measures the ini- 
tial number of phonons. Vice versa, if the component 2 
is not empty initially, the Raman beams transfer atoms 
(Cfe+K.) from the state 2 to the level 1 with simultaneous 
emission of an equal number of phonons (a^). 

The energy-momentum balance (|14|) is a bit more 
complicated than in the previous case without interac- 
tions, but exhibits a similar direction-degeneracy, which 
can again be eliminated by considering effectively one- 
dimensional condensates. In the phonon limit (A ^ ^ 
and (5 <C /i), we obtain a unique solution for the phonon 
energy LOk ^ 5 which allows us to address single modes 
with suitably tuned Lasers. If we choose 5 and n to lie 
on the phonon dispersion curve 5 = uj{k), we annihilate 
one phonon with energy 6 and momentum k and create 
one particle in the component r = 2 in the ground state. 
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With sufficiently long pulses leading to a good energy 
resolution, it should be possible to "see" the discrete na- 
ture of the phonon spectrum, i.e., to address single (or a 
few) phonon modes. In order to annihilate all phonons 
in the r = 1 condensate with a given energy/momentum 
and to transfer the same number of atoms to the r — 2 
component, we apply an effective Raman 7r-pulse with 
the duration [cf. Eq. l(TK|l ] 
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where the « sign applies to the phonon limit. 

Of course, the approximations used in the presented 
derivations must be checked for a potentially realistic 
set of experimental parameters. Let us assume a speed 
of sound of a few millimeters per second and a healing 
length around one micrometer. In this case, the wave- 
length A of the phonons to be detected would typically 
be several micrometers 1/k = 0(10 /xm) and their fre- 
quency a few hundred Hertz 5 = 0(100 Hz). 

Using Lasers in the optical range 0(10^^ Hz), the 
large detuning A depends on the atomic level struc- 
ture and would be a little bit below this value, say 
A = 0(10^3 _ 1014 Hz). With quite moderate Rabi fre- 
quencies — 0(10^* — 10'' Hz), we can achieve an effec- 
tive Raman transition rate \/ [ijb 17^/ A of a few tens of 
Hertz. Consequently, the duration of the effective Raman 
TT-pulse in Eq. (|16|l would be of the order of hundred mil- 
liseconds T — 0(100 ms) leading to a energy resolution 
of circa ten Hertz. In view of the aforementioned pa- 
rameters, the assumptions and approximations (e.g., the 
adiabaticity A ^ fi) used in the derivation are reason- 
ably well justified. The major constraint is given by the 
energy resolution of the effective Raman 7r-pulse peaked 
around 8 ^ 0(100 Hz) ± 0(10 Hz). Apart from a few 
excitations (i.e., phonons), the beams illuminate many 
atoms in the ground state (zero energy) and one has to 
make sure that the probability of transferring an atom 
from the ground state of the condensate in component 
r = 1 into the state r = 2 is small enough. Thus the 
negative-frequency tail of the Fourier transform of the 
pulse (which is peaked around b in frequency space) must 
be suppressed accordingly 

With the ability of measuring a few low-energy 
phonons, it might become possible to observe some of 
the exotic quantum effects mentioned in the Introduc- 
tion. The analogue of cosmological particle production 
is probably easier to realize experimentally than Hawk- 
ing radiation since it can be done with a condensate at 
rest and a practically unlimited measurement time (after 
varying Cs via a Feshbach resonance, for example). The 
same advantage applies to a small wiggling stirrer in the 
condensate, which would act as a point-like non-inertial 
scatterer and generate the analogue of moving-mirror ra- 
diation [^, which can be interpreted as a signature of 



the Unruh effect. In contrast, the detection of the Hawk- 
ing radiation requires either a flowing condensate or a 
motion of the horizon via a space-time dependent sound 
velocity Cg (t, x) , cf. |l4l| . Apart from measuring this strik- 
ing effect, these experiments may also shed light onto 
the trans-Planckian problem, i.e., impact of the short- 
range physics on the long-wavelength Hawking radiation: 
Even though the Hawking effect seems to be quite robust 
against modifications of the dispersion relation at short 
wavelengths (such as the Bogoliubov dispersion in Bose- 
Einstein condensates) only very little is known about the 
impact of interactions, see, e.g., [l5| . 
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